雖然特徵多項式 $p(\lambda) = \det(A - \lambda I)$ 是定義特徵值的理論基礎，但對於高維系統而言，它在數值上「條件不良」且計算效率低下。在實際應用中——例如求解波傳播的 Sturm-Liouville 系統——多項式根對係數擾動極其敏感，因此直接展開成為次優選擇。
從連續波到離散矩陣
弦或膜的振動由波動方程支配：
$$\rho(x) \frac{\partial^2 v}{\partial t^2}(x, t) = \frac{\partial}{\partial x} \left[ p(x) \frac{\partial v}{\partial x}(x, t) \right]$$
為求得解 $v(x, t) = \sum_{k=0}^{\infty} c_k u_k(x) \cos \sqrt{\lambda_k}(t - t_0)$，我們必須求解 Sturm-Liouville 系統：
$$\frac{d}{dx} \left[ p(x) \frac{du_k}{dx}(x) \right] + \lambda_k \rho(x) u_k(x) = 0$$
離散化複雜性
將算子離散化會產生如 $Aw = -0.04 \frac{\rho}{p} \lambda w$ 之類的矩陣方程。對於 $4 \times 4$ 的三對角矩陣，$p(\lambda)$ 尚可處理。然而，隨著網格細化（$n$ 增加），我們將面臨兩大障礙：
- 阿貝爾-魯菲尼限制： 當 $n \ge 5$ 時，多項式的根不存在代數解。
- 捨入敏感性： 在高維系統中，一個元素的小數點後第 $10^{-10}$ 位的微小變化，可能導致特徵值移動數個數量級（威爾金森多項式現象）。
數值必要性與專業函式庫
專業的數值函式庫（IMSL、NAG）避免使用原始特徵多項式，而是採用迭代方法進行近似：
- IMSL 函式庫： 使用線性最小二乘法、三次樣條插值及快速傅里葉變換。
- NAG 函式庫： 採用最小二乘多項式逼近與 $l_1/l_{\infty}$ 范數擬合。
在近似系統 $\lambda_i = 1 + 4\alpha\left(\sin \frac{\pi i}{2m}\right)^2$ 的特徵值時，我們依賴離散最小二乘法與迭代發現，而非根尋找。
🎯 理論工具與數值危機
特徵多項式 $p(\lambda) = \det(A - \lambda I)$ 對於證明至關重要，但在計算中卻具有危險性。物理學中的實務特徵值問題透過迭代變換（如 QR 法）解決，以確保穩定性。